Quantum simulation of classical thermal states 
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We establish a connection between ground states of local quantum Hamiltonians and thermal 
states of classical spin systems. For any discrete classical statistical mechanical model in any spatial 
dimension, we find an associated quantum state such that the reduced density operator behaves 
as the thermal state of the classical system. We show that all these quantum states are unique 
ground states of a universal 5-body local quantum Hamiltonian acting on a (polynomially enlarged) 
system of qubits arranged on a 2D lattice. The only free parameters of the quantum Hamiltonian 
are coupling strengthes of two-body interactions, which allow one to choose the type and dimension 
of the classical model as well as the interaction strength and temperature. 
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Introduction. — Classical statistical mechanical models 
are widely studied in physics and have found applications in 
various contexts, ranging from magnetism and neural net- 
works [1J to lattice gauge theories describing fundamental 
interactions in nature [2]. Thereby models with arbitrary 
underlying geometry and different symmetry are included, 
e.g. Ising or Potts models [3] that have a global symmetry 
or abelian discrete lattice gauge theories with local symme- 
tries. Many of these models show rich phase diagrams and 
are known to be related to mathematical problems with 
large computational complexity [4 . Their inherent com- 
plexity and limited experimental accessability, e.g. due to 
a spatial dimension larger than three, hinders in some cases 
our understanding of these systems. 

Here we propose a quantum simulator for thermal states 
of all discrete classical spin systems in arbitrary spatial di- 
mension. We show that there exists a universal, local 2D 
quantum Hamiltonian whose unique ground state contains 
the thermal state of all discrete classical spin systems, in- 
cluding systems with long-ranged and many-body interac- 
tions in arbitrary geometries. The choice between different 
classical models and parameters is provided by the coupling 
strength of two-body interactions in the 5-body quantum 
Hamiltonian. For all reasonable discrete classical models, 
the size of the quantum system is only polynomially en- 
larged. This opens the possibility to study thermal states 
of classical spin models in arbitrary dimension using a well- 
controlled 2D quantum system. 

The result is inspired by recently established connections 
between classical spin systems and quantum systems (5j-[9] , 
which have been utilized e.g. to identify complete classical 
spin models such as a 4D lattice gauge theory, i.e. mod- 
els from which the partition function of all other classical 
models can be obtained |MT0|. Here we prove a similar 
completeness result, where we however provide a link be- 
tween all discrete classical spin models and ground states of 
local quantum Hamiltonians. The relation is not only at the 
level of the partition function, but the thermal state of the 
classical spin system is indeed contained in the quantum 
system in a physical sense. 

To show this relation, we proceed in three steps: (i) We 
first introduce a new mapping between classical spin sys- 



tems and quantum states, where for each discrete classical 
spin model a (polynomially enlarged) quantum state is de- 
fined that contains the thermal state of the classical model 
in a certain subsystem, i.e. by simply tracing out or ignor- 
ing the additional quantum particles, (ii) We then make 
use of results from measurement-based quantum computa- 
tion |11[I12|. which shows how to obtain arbitrary quantum 
states from a universal state, the so-called 2D cluster state 
|13j . We use these insights to write the quantum state cor- 
responding to any classical spin system as a polynomially 
enlarged deformed 2D cluster state, where (non-unitary) 
deformations act on individual qubits. (iii) Finally we show 
that these deformed cluster states can always be written as 
the ground state of a universal, local 5-body 2D quantum 
Hamiltonian, where only coupling strengths of two-body 
interactions appear as free parameters. The latter tech- 
nique which allows one to (approximately) write all effi- 
ciently preparable quantum states as unique ground states 
of local Hamiltonians might also find applications in other 
contexts. 

Mappings between classical spin systems and quantum 
states. — Consider a classical spin model with N 2-level 
systems s — (si, S2, ■ ■ ■ sjv), *t € {0, 1} which contains all 
many-body interactions of Ising type between groups of at 
most d spins. For notational simplicity we will here restrict 
attention to d = 3 but the generalization to higher d is 
straightforward. The system is described by the Hamilto- 
nian 

H(s) = 2J Ja(s a ) + ^ J a.b{Sa, Sb) + ^ J abc(s a , S b , S c ). 
a ab abc 

The function J a specifies an arbitrary local energy term 
(magnetic field) depending on s a and J a b, J a b c associate to 
each spin configuration of a subset of 2 or 3 spins an en- 
ergy that only depends on the parity of the involved spins, 
e.g. J ab {s a ,s b ) = (-l) Sa ® Sb J a b (which we call an Ising-type 
interaction). The Hamiltonian H contains all possible 1-, 
2- and 3-body interactions of this kind. Hence, by setting 
suitable couplings to zero, Ising-type models on arbitrary 
lattices with 3-body interactions can be obtained. 

Interestingly, Ising-type Hamiltonians H are in fact suf- 
ficient to treat also non-Ising models with arbitrary /c-body 
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interactions and g-level systems i.e. H defines a complete 
classical spin model, as shown in Roughly speaking, 
one notes that q-level spins can be encoded into |~log 2 q] 
2-level spins, and that arbitrary fc-body interactions of q- 
level spins can be decomposed into (all possible) m-body 
Ising-type interactions with m < kx \log 2 q~\ ■ The required 
Ising interaction strengths can be determined by solving 
a system of linear equations [9]. As long as q and k are 
bounded, this type of encoding is efficient as the overhead 
in additional qubits scales polynomially with N . 

We now define a quantum state that "contains" the 
thermal state of the classical system H . We proceed in 
two steps. First we associate a qubit with every classi- 
cal spin (vertex particles) and with every interaction term 
Jab(s a , Sb), J a bc{s a , Sfej s c ) (interaction particles) and we in- 
troduce a so-called d-clique state that encodes the interac- 
tion pattern (see also [5]) of H: 

\<p) = J2 < 8> i Sa ® Sb ® Sc > ® i s ° ® sfc > ® I Sq >- ( 2 ) 

s abc ab a 

Notice that \(p) is a stabilizer state. In the second step, we 
diverge from the approach of |8J and add information about 
interaction strengths and inverse temperature /3 = (fc^T) -1 
via (non-unitary) deformations. We define diagonal matri- 
ces A a , A ab , Aabc, where e.g. A ah = diag(e -/3J ° b/2 , e^ Jab/2 ), 
and analogous for A a ,A abc . We use the short-hand nota- 
tion A = o6c A abc afc A ab (g) Q A a for the tensor product 
of all single-qubit deformations. We consider the deformed 
stabilizer state 

\<p A ) - A\<p)/Vz. (3) 

Here the partition function Z — ^2 s e~^ H ^ appears as a 
normalization. What is more, if one traces out all inter- 
action particles, the resulting density matrix of the vertex 
particles has as diagonal entries the Boltzmann weights of 
H [H], p = Tr abtabc I^aX^aI with 

(s\p\s) = e-" H ^/Z. (4) 

Therefore all classical quantities such as spin correlations 
or energies are contained in p and can be read out using 
diagonal quantum observables. 

We remark that the states \ip\) can be written as ground 
states of certain d-body quantum Hamiltonians. However, 
these Hamiltonians may be highly non-local and may con- 
tain <3(iV)-body interactions terms. In the following we will 
show how to relate this state to a local two-dimensional 
Hamiltonian. To this aim it is crucial that the deforma- 
tions in |<^a) are single-qubit operators, which is ensured 
by the specific mapping above making use of only Ising type 
interactions. 

Universal quantum state. — We now show that, inde- 
pendent of the choice of parameters, any state can 
be written to arbitrary accuracy as a deformed 2D-cluster 
state with invertible single-qubit deformation operators. 
The 2D cluster state \C) [13J is a universal resource for 
measurement-based quantum computation (MQC), in the 
sense that any quantum state can be obtained from a suffi- 
ciently large 2D-cluster state by means of single-qubit mea- 
surements [TTI [H] . The 2D-cluster state is associated with 



a 2D rectangular lattice by placing qubits on the vertices, 
and is defined as 

\C)=U\+)® M , U= H U ab (5) 

(a,b)eE 

where |+) = ^(|0) + |1», U ab = diag(l, 1, 1, -1) and E 

are the edges of the 2D lattice. The universality of the 2D 
cluster state implies that there exists single-qubit states 
\aj a ) such that 

(®|w«XWa|®/)|C) = ^ r (g)|w„)®b), (6) 
\a£A ) 2 2 aeA 

i.e. by projecting out all particles in A, the stabilizer state 
1^) can be generated at the remaining particles. Notice 
that this method is constructive and efficient, i.e. a poly- 
nomially enlarged 2D cluster state suffices and only Pauli 
measurements are involved since \ip) is a stabilizer state. In 
order to relate the resulting quantum state to a local Hamil- 
tonian, we now replace the projection operators |u; a )(aj a | 
by invertible operators Q a = (1 — e) |cj a )(w a | + e 
where a polynomially (in N) small e suffices to guarantee 
an approximation of the desired state \<p) with polynomial 
accuracy |15j . To be precise, the resulting state after ap- 
plying imperfect projectors Q a at each of the \A\ qubits is 
given by a sum of 2' A ' terms, corresponding to all possible 
combinations where at n positions the desired projection 
|w a )(u; a | has been applied, while at |^4| — n positions the 
(wrong) projector Iw^Xw^l was applied. Every term has a 
weight (1 — e } n e \ A \— n i while each of the resulting states has 
the same norm, independent of which projection operators 
have been applied. The latter point is a special feature of 
2D cluster states that also occurs in measurement-based 
quantum computation, where outcomes of projective mea- 
surements in any basis occur with equal probability at any 
stage of the computation. The reason is that reduced den- 
sity operator of each qubit is given by the identity due to 
the entanglement with the remaining (unmeasured) qubits. 
We can then lower bound the fidelity of the resulting state 
as compared to the one obtained by the perfect projections 
by F > (1 — e)' A ', where we assume that only the term 
where all projections lead the desired outcome contribute 
to the fidelity. Using the notation Q = ^> aeA Q a , it follows 
that the deformed stabilizer state |(/?a) can be written to 
arbitrary accuracy as a deformed 2D cluster state 

|C n , A > =M® A|C) w ® |w„) ® |^a>, (7) 

aeA 

where we have omitted the normalizations. Notice that 
there are 3 different groups of qubits: group A are auxiliary 
qubits, which are deformed by to form the state |<^a) on 
the remaining qubits. The second group B are all interac- 
tion qubits corresponding to Ising-type interactions, where 
the deformation is determined by the interaction strengths 
J ab and J abc and the inverse temperature /?. Finally, the 
group of vertex qubits C is acted upon by deformations 
determined by the local fields J a and by (3. At the same 
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time, C is the group on which the thermal state of the 
corresponding classical system in generated. 

Ground states of local quantum Hamiltonians. — We now 
show that the deformed cluster state ICq.a) can always be 
written as the unique ground state of a universal local 5— 
body Hamiltonian of a 2D quantum system. 

The deformations f2 and A are treated in different ways. 
We first consider the partially deformed state \C\) = 
I ® A\C). Note that A (which acts on groups B and C) 
is a diagonal matrix which hence commutes with the phase 
gates U so that |C A > = U[I ® A\\+)® N . Remark that 
for every i £ A the state is the unique ground state 
of —Xi with ground state energy — 1. Moreover each of 
the local deformations transforms \+)k with k G B U C in 
e -/3J / 2 |0) + e? J l 2 \\). This state is the unique ground state 
of a single-qubit Hamiltonian — X k — ^ k Z for some suitable 
7^ with ground state energy, say, E k . For every a consider 
the 5-local operator K a — X a Y[ Zt where the product is 
over all neigbors b of a in the 2D lattice; the K a are the 
standard stabilizer operators of the cluster state. Using 
that UX a W = K a and UZ a U^ = Z a , we find that |C A > is 
the unique state satisfying 

{I-Ki)\C A )=0 {-K k - lk Z k -E k I)\C A )={) (8) 

for every i e A and fc e BUC, and is hence the zero-energy 
ground state of the Hamiltonian containing the sum of all 
these terms. 

We now deal with the deformations f2 acting on group 
A. The operator Jl is in general not diagonal, so a differ- 
ent treatment is necessary. First, let f2* denote the tensor 
products of all fi which act on i and its nearest-neighbors 
in the 2D lattice. Since \C\) is a zero energy ground state of 
I — Ki, it follows that |Co,a) is a zero energy ground state 
of Pi = [n*]-t(/ - K^fi 1 }- 1 for every i € A [16]. This 
particular choice of transformation guarantees that the re- 
sulting Hamiltonian Pi is hermitian, and remains 5-body 
|17j . Notice that here it is important that the inverse of 
£l l exists, which is the reason why we approximate pro- 
jection operators by invertible operators flf, in Eq. |7]). 
Second, we consider the terms H k = —K k — ^ k Z k — E k I, 
k G B U C of the Hamiltonian, where in each term only one 
particle n k £ A is affected by a deformation in group A 
(see Fig. lb). Hence |Cn,A) is a zero-energy ground state 
of Q k = [ft W) J _ t.Hfc[ft n) J _1 . Putting everything together, 
we find that |Cn a) is the unique ground state of the 2D 
5-body Hamiltonian 

n = /2 p *+ E o*- ( 9 ) 

ieA keBuc 

Notice that term —j k Z k in H k is transformed into a two- 
body interaction 7fc([f2n fc ] [finj -1 ) Z k . The strength 
of this fixed two-body term is determined by the parame- 
ters of the deformation A (and hence the parameters of the 
corresponding classical model) . By the completeness of the 
classical model , one can thus change between all classi- 
cal models (which can correspond to arbitrary geometries 
and also to local or global symmetries) by simply adjusting 
the strength of these two-body interactions. Notice that in 



this way we obtain a constructive method to realize thermal 
states of arbitrary classical models in any spatial dimension 
as ground states of a 2D local quantum Hamiltonian. This 
proves the main result of this paper. 

Note that a reduction to a local 4-body Hamiltonian is 
possible by using a universal graph state corresponding to 
a 2D hexagonal lattice instead of the 2D cluster state. 

Direct constructions. — It is straightforward to repeat the 
above argument starting from Ising-type models on arbi- 
trary lattices rather then the complete model 0. This 
involves a stabilizer state similar to \(p) which now encodes 
the interaction pattern; see Fig [T] for an example. The 
application of the deformations A is analogous, as is the 
mapping to a deformed cluster state which is then found to 
be a ground state of a 2D 5-body Hamiltonian as well. 




FIG. 1: (a) Illustration of mapping for a classical spin system 
on a triangular lattice with nearest neighbor 2-body and 3-body 
plaquette interactions. Quantum particles are associated with 
vertices (green, large), nearest neighbor pairwise interactions 
(red, medium) and 3-body plaquette interactions (blue, small). 
The classical thermal state is generated at the vertex particles 
(green, large), (b) 2D Cluster state with groups A, B, C. De- 
formations on group A allow one to generate the state \<p) on 
system B U C (Fig. 1(a)), while deformations on group B and 
C allow one to chose interaction strengthes and temperature. 

Discussion and Conclusion. — We have shown that there 
exists a universal family of 2D 5-body quantum Hamil- 
tonians Eq. (|9| whose unique ground states include the 
thermal states of all classical spin models in arbitrary spa- 
tial dimension. For all classical spin models with g-level 
spins and fc-body interactions where q and k are bounded 
-including edge models, vertex models as well as lattice 
gauge theories-, the size of the corresponding quantum 
system is polynomially enlarged, and the classical thermal 
state is generated at a subsystem of the quantum spins. 
The parameters and type of the classical model are deter- 
mined by the strength of fixed two-body interactions in the 
quantum Hamiltonian. 

Apart from this general result, two techniques applied 
in this proof are of particular relevance and might find ap- 
plications elsewhere. On the one hand, the new mapping 
relating classical models with deformed stabilizer states al- 
lows one to study the feature of classical models by investi- 
gating the (entanglement) properties of the corresponding 
deformed stabilizer states. On the other hand, the tech- 
nique of relating a deformed cluster state to a ground state 
of a local Hamiltonian also allows one to write quantum 
states as ground states of a local quantum Hamiltonian. It 
would be interesting to see if the possibility to efficiently 
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prepare these states via MQC is reflected in the properties 
of the corresponding Hamiltonian, e.g. its gap. 

We also point out that our results lead to a quantum 
simulator for classical thermal states. This is either based 
on the direct generation of the ground states of the local 
quantum Hamiltonian Eq. (J9j> , e.g. by a cooling process or 
by adiabatic quantum computation, or by the generation of 
the states \ipj^) by some other means, e.g. using a quantum 
circuit. For many classical models this yields an efficient 
quantum algorithm, e.g. for (inhomogenous) classical mod- 
els with pairwise interaction pattern given by a graph with 
bounded rank width, e.g. ID systems or tree graphs where 
the corresponding quantum state is a tree tensor network 
state [20J, or the the 2D Ising model without fields (see 
[7] for more examples). Notice however that one cannot 
expect that an efficient preparation of the corresponding 
quantum state is possible for all models, as this would im- 
ply the solution of NP-hard problems using our method. 
One can rather expect that in such cases the quantum state 
preparation becomes inefficient. It would be interesting to 
study the connection between computational complexity of 
the classical spin models and the possibility to efficiently 
generate the corresponding deformed cluster states, in par- 
ticular the conditions under which an efficient generation 
of the quantum state is possible, e.g. by using recently de- 
veloped algorithms to efficiently prepare certain projected 
entangled pair states [21]. The gap in the Hamiltonian Eq. 
([9]) plays a key role in this context, and will be subject of 
further study. 



We finally remark that a similar result to ours might 
be obtained making use of connections between classical 
thermal states and quantum states pointed out in [S] , each 
of which can be written as a ground state of a (non- 
local) quantum Hamiltonian. A reduction to 2D two-body 
Heisenberg Hamiltionian with only local fields as control 
parameters could be achieved by making use of gadgets 
constructions of Refs. [T51 [TU] . This method however has 
the disadvantage that already for models with long-ranged 
two-body interactions, i.e. where each spin interacts with 
O(N) other spins, the corresponding quantum Hamiltonian 
is iV-body. This leads to an exponential overhead in the 
number of auxiliary particles in the gadget construction. 
In addition, the gadget constructions put stringent require- 
ments on the necessary accuracy of the control parameters 
of the Hamiltonians |181 119) . Notice that gadget construc- 
tions can also be applied as a final step in our construction, 
reducing the universal quantum Hamiltonian Eq. Q to an 
effective two-body Hamiltonian. In addition, the two-body 
interactions controlling the type and parameters of the clas- 
sical model can be reduced to a local single-particle mag- 
netic field together with fixed two-body interactions, i.e. 
the change between different classical models is achieved 
by varying local fields only. 
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